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1.  Introduction 


Missile  concepts  with  forward  control  fins,  or  canards,  have  been  used  for  many  years. 

However,  previous  studies  have  shown  that  concepts  with  canards  can  suffer  from  adverse 
induced  rolling  moments  [1-4].  The  use  of  grid  fins,  or  “lattice  controls,”  for  the  tail  control 
surfaces  instead  of  conventional  planar  fins  was  recently  proposed  as  a  possible  remedy  for  the 
roll  control  problems  [5]. 

A  grid  fin,  shown  in  section  2,  is  an  unconventional  lifting  and  control  surface  that  consists  of  an 
outer  frame  supporting  an  inner  grid  of  intersecting  planar  surfaces  of  small  chord.  The 
aerodynamics  of  grid  fins  has  been  investigated  since  1985  by  the  U.S.  Army  Aviation  and 
Missile  Command  (AMCOM),  Redstone  Arsenal,  AL  [6-8].  Studies  have  shown  that  when 
compared  to  conventional  planar  fins,  grid  fins  have  advantages  such  as  effective  aerodynamic 
control  at  high  a  and  high  Mach  number,  a  small  hinge  moment,  compact  storage,  improved  yaw 
stability,  attenuated  body-vortex  interference,  and  improved  roll  control  of  canard-controlled 
munitions  [5-7, 9].  The  primary  disadvantage  of  the  grid  fin  concept  is  a  higher  drag  than 
conventional  planar  fins,  which  can  be  minimized  with  the  proper  design  [8].  Computational 
fluid  dynamic  (CFD)  techniques  to  calculate  the  viscous  flow  around  grid  fins  were  recently 
demonstrated  [10, 11]. 

In  this  study,  the  CFD  techniques  demonstrated  previously  [10, 11]  were  applied  to  study  the  roll 
control  effectiveness  of  both  planar  and  grid  fins  on  a  generic  canard-controlled  missile  shape. 
The  numerical  calculations  were  validated  with  wind  tunnel  data  from  an  experimental  study 
from  AMCOM  and  the  former  Defence  Research  Establishment  Valcartier  (DREV),  now  the 
Defence  Research  and  Development  Canada-Valcartier  (DRDC-V)  [5].  The  objectives  of  the 
CFD  study  were  to  provide  insight  into  the  fluid  dynamic  phenomena  that  cause  roll  reversal  on 
canard-controlled  missiles  with  planar  tail  fins  and  permit  the  retention  of  roll  control  when  grid 
tail  fins  are  used  on  canard-controlled  missiles.  The  CFD  simulations  provide  visualizations  of 
the  flowfield  that  aid  the  understanding  of  the  flow  physics  responsible  for  the  adverse  forces 
and  moments  observed  in  the  wind  tunnel  and  calculated  results. 


2.  Computational  Approach 


2.1  Geometry  and  Simulation  Parameters 

The  investigation  used  CFD  to  determine  the  flowfield  and  aerodynamic  coefficients  on  a 
16-cal.,  four-finned,  generic  canard-controlled  missile.  This  study  followed  the  experimental 
wind  tunnel  investigation  performed  at  DREV  [5].  The  DREV  wind  tunnel  is  an  intermittent,  in¬ 
draft  wind  tunnel  with  a  0.6  x  0.6-m  test  section.  In  this  type  of  tunnel,  the  air  flows  from  an 
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atmospheric  pressure  tank  to  a  vacuum  tank,  and  the  Reynolds  number  is  lower  than  free-flight 
values  at  high  Mach  numbers.  The  wind  tunnel  Reynolds  number  ranges  from  about 
1.56  x  107  m'1  at M  =  1.15  to  4.7  x  106  m'1  at  A/=  4.0. 

The  wind  tunnel  model  geometry  was  used  in  the  CFD  study.  Four  canards  on  the  ogive  were  in 
line  with  the  fins.  Two  fin  types  were  investigated:  conventional  planar  fins  and  grid  fins. 

Figure  1  and  Figure  2  show  the  geometry  for  the  planar  fin  and  grid  fin  cases,  respectively.  The 
missile  was  16  cal.  long  with  a  hemispherical  nose  on  a  3.7-cal.  truncated  tangent  ogive.  The 
canard  midchord  was  located  0.96  cal.  from  the  missile  nose,  and  the  fin  midchord  was  located 
1 .5  cal.  from  the  missile  base.  The  canards  (Figure  3a)  had  a  double-wedge,  trapezoidal 
planform  with  a  span  of  0.37  cal.,  a  root  chord  of  0.36  cal.,  a  tip  chord  of  0.13  cal.,  a  midchord 
root  thickness  of  0.03  cal.,  and  a  taper  ratio  of  1.48.  The  planar  fins  had  a  double-wedge, 
rectangular  planform  with  a  span  of  0.78  cal.,  a  chord  of  0.65  cal.,  and  a  midchord  thickness  of 
0.03  cal.  The  grid  fin  structure  (Figure  3b)  consisted  of  23  cubic  and  12  prismatic  cells  with  a 
span  of  0.74  cal.,  a  chord  of  0.10  cal.,  and  a  thickness  of  0.46  cal.  The  web  thickness  between 
the  grid  fin  cells  was  0.003  cal.  The  canard  and  fin  characteristics  are  summarized  in  Table  1 . 


Figure  3.  (a)  Missile  nose  with  canards  at  8  =  0°  and  (b)  front  view  of  missile. 


Table  1.  Canard  and  fin  characteristics. 


Control  Type 

Span 

(cal.) 

Root  Chord 
(cal.) 

Tip  Chord 
(cal.) 

Root  Thickness 
(cal.) 

Taper  Ratio 

Web 

(cal.) 

Canard 

0.37 

0.36 

0.13 

0.03 

1.48 

- 

Planar  fin 

0.78 

0.65 

0.65 

0.03 

1.0 

— 

Grid  fin 

0.74 

0.10 

0.10 

0.46 

1.0 

0.003 

The  analyses  were  performed  at  two  Mach  numbers:  M  =  1.5  and  3.0;  two  canard  deflections: 

8  =  0°  and  10°;  and  six  angles  of  attack:  a  =  0°,  2°,  4°,  6°,  8°,  and  10°.  The  DREV  wind  tunnel 
conditions  were  used  in  this  study.  For  M  =  1.5,  the  freestream  conditions  were  a  Reynolds 


number  of  1.4  x  107  m"1,  a  static  temperature  of  206  K,  and  a  static  pressure  of  2.6  x  104  Pa.  For 
M=  3,  the  ffeestream  Reynolds  number  was  7.6  x  106  m"1,  the  static  temperature  was  107  K,  and 
the  static  pressure  was  2.77  x  103  Pa.  The  model  reference  diameter  was  30  mm,  and  the 
moment  reference  point  was  10.63  cal.  aft  of  the  missile  nose.  The  simulations  were  performed 
with  the  missile  in  the  cruciform  (+)  configuration.  The  DREV  wind  tunnel  data  ranged  from 
-4°  <  a  <  +15°.  In  the  8=10°  case,  all  four  canards  were  deflected  in  the  same  direction, 
intended  to  give  a  positive  roll,  which,  by  convention,  was  clockwise  when  viewed  from  the  rear 
of  the  missile. 

2.2  Solver 

Steady-state  calculations  were  used  to  compute  the  flowfield  using  the  commercial  CFD  code, 
FLUENT  version  5.5  [12].  The  implicit,  compressible,  unstructured-mesh  solver  was  used.  The 
three-dimensional  (3-D),  time-dependent,  Reynolds-Averaged  Navier-Stokes  (RANS)  equations 
are  solved  using  the  finite  volume  method: 

—  |W4F  +  j[F-G]  dA  =  JH4F,  (1) 

at  v  y 

where 


p 

pv 

0  ' 

pu 

pvu  +  pi 

pv 

.,  F  =  < 

pvv  +  p] 

>,  G  =  < 

U/ 

pw 

pvw+ pk 

u, 

PE. 

pvE  +  pv 

Tuvj  +q 

The  inviscid  flux  vector  F  is  evaluated  by  a  standard  upwind  flux-difference  splitting.  In  the 
implicit  solver,  each  equation  in  the  coupled  set  of  governing  equations  is  linearized  implicitly 
with  respect  to  all  dependent  variables  in  the  set,  resulting  in  a  block  system  of  equations.  A 
block  Gauss-Seidel,  point  implicit  linear  equation  solver  is  used  with  an  algebraic  multigrid 
method  to  solve  the  resultant  block  system  of  equations.  The  coupled  set  of  governing  equations 
is  discretized  in  time  and  time  marching  proceeds  until  a  steady-state  solution  is  reached.  In  the 
implicit  scheme,  which  was  used  in  this  study,  an  Euler  implicit  discretization  in  time  is 
combined  with  a  Newton-type  linearization  of  the  fluxes. 

A  modified  form  of  the  k-s  two-equation  turbulence  model  was  used  in  this  study.  Called  the 
“realizable”  k-s  model  in  FLUENT,  it  differs  from  the  standard  k-s  model  in  that  it  contains  a 
new  formulation  for  the  turbulent  viscosity  and  a  new  transport  equation  for  the  dissipation  rate, 
which  was  derived  from  an  exact  equation  for  the  transport  of  the  mean-square  vorticity 
fluctuation  [12].  The  term  “realizable”  means  that  the  model  satisfies  certain  mathematical 
constraints  on  the  Reynolds  stresses  consistent  with  turbulent  flow  physics.  The  realizable 
k-s  model  has  shown  substantial  improvements  over  the  standard  k-s  model  where  flow  features 
include  strong  streamline  curvature,  vortices,  and  rotation  [12]. 
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2.3  Computational  Mesh  and  Boundary  Conditions 

The  geometry  and  unstructured  mesh  were  generated  using  the  preprocessor,  GAMBIT,  supplied 
in  the  FLUENT  software  suite.  Canard  deflection  and  angle  of  attack  precluded  the  use  of 
symmetry  or  periodicity,  so  a  full  3-D  mesh  was  required.  In  generating  the  meshes,  boundary 
layer  mesh  spacing  was  used  near  the  missile  body  and  fin  surfaces.  The  two-layer  zonal  model 
was  used  for  the  near- wall  equations  and  the  first  point  off  the  surface  (cell  center)  was  about 
7.0  x  10~5  cal.  from  the  surface,  chosen  to  give  a  y+  value  of  about  1 .0.  All  mesh  stretching  was 
kept  at  or  below  1 .25.  About  144  cells  were  used  on  the  missile  body  in  the  circumferential 
direction,  with  this  value  increased  in  the  grid  fin  region. 

An  all-hexahedral  (hex)  mesh  was  used  for  the  planar  fin  case,  with  the  total  number  of  cells 
about  4.9  M.  O-grid  type  meshes  were  generated  around  the  canards  and  planar  fins.  Figure  4a 
and  Figure  4b  show  the  surface  mesh  in  the  tail  region  and  ogive  region,  respectively.  A  hybrid 
hex  and  tetrahedral  (tet)  mesh  was  used  for  the  grid  fin  case,  with  the  total  number  of  cells  about 
15.7  M.  The  first  13  cal.  of  the  missile  was  meshed  with  the  same  hex  mesh  as  in  the  planar  fin 
case.  A  tet  mesh  was  used  in  the  tail  region  (Figure  5a)  in  order  to  mesh  the  complicated  grid  fin 
structure.  Layers  of  triangular  prisms  were  used  on  the  body  to  capture  the  boundary  layer 
(Figure  5b).  The  tet  mesh  in  the  tail  is  matched  to  the  forward  hex  mesh  via  a  layer  of  pyramidal 
cells.  Due  to  meshing  constraints,  prism  layers  were  not  used  around  the  grid  fins,  so  the  spacing 
of  the  first  point  off  the  grid  fin  surfaces  was  larger  than  optimum.  Postprocessing  of  the  runs 
showed  that  they1"  values  on  the  missile  body,  canards,  and  planar  fins  were  less  than  1.0.  They+ 
value  on  the  grid  fins  ranged  from  about  2-15,  with  most  of  the  grid  fin  surface  near  the  higher 
value.  Although  these  values  are  not  optimum  for  calculating  the  boundary  layer  properties  on 
the  grid  fins,  there  was  little  apparent  effect  on  the  missile  aerodynamic  coefficients. 

The  base  flow  was  not  simulated  in  these  calculations,  so  the  mesh  ended  at  the  base  of  the 
missile.  For  supersonic  flow,  this  has  been  found  to  be  reasonable  methodology  if  any  tail  fin 
surfaces  are  not  located  too  close  to  the  base  of  the  missile.  The  computational  domain  extended 
about  4  cal.  from  the  missile  body.  An  outflow  boundary  condition  was  used  downstream,  a 
pressure  inflow  (with  freestream  conditions)  boundary  condition  was  used  upstream,  and  a  far- 
field  pressure  (nonreflecting)  boundary  condition  was  used  for  the  outer  boundary.  A  no-slip 
wall  boundary  condition  was  used  for  all  solid  surfaces. 

2.4  Solution  Methodology 

The  simulations  were  performed  in  parallel  on  Silicon  Graphics,  Inc.  (SGI)  Origin  2000  and 
3800  machines.  The  simulations  were  run  with  a  maximum  Courant-Friedrich-Lewy  (CFL) 
number  of  10  for  the  planar  fin  cases  and  5  for  the  grid  fin  cases.  Each  case  was  started  with  a 
lower  CFL  value  (usually  about  1.0)  and  ramped  up  to  the  maximum  during  the  first  few 
hundred  iterations  of  the  simulation.  The  calculations  took  -300-600  ps/cell/iteration,  using 
28  processors  for  the  planar  case  and  64  or  96  processors  for  the  grid  fin  case.  Convergence 
was  determined  by  tracking  the  change  in  the  flow  residuals  and  the  aerodynamic 
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Figure  4.  Views  of  computational  mesh:  (a)  tail  region  of 
planar  fin  case  and  (b)  ogive  region. 

coefficients  during  the  solution.  The  solution  was  deemed  converged  when  the  flow  residuals 
had  stabilized  and  the  aerodynamic  coefficients  were  changing  less  than  0.5%  after  the  last  100 
iterations.  The  aerodynamic  coefficients  converged  in  about  1000-1500  iterations. 

The  planar  fin  cases  were  run  with  the  double  precision  solver.  To  reduce  the  computer  memory 
required,  the  grid  fin  cases  were  run  with  the  single  precision  solver.  A  planar  fin  case  (a  =  4°, 
8=10°)  was  repeated  with  the  single  precision  solver  to  determine  if  there  was  an  effect  on  the 
calculated  solution.  The  six  aerodynamic  coefficients  calculated  with  the  single  precision  solver 
were  within  0.5%  of  those  calculated  with  the  double  precision  solver.  Four  of  the  coefficients 
were  within  0.1%.  Therefore,  the  single  precision  solver  was  adequate  for  these  computations. 


(a) 


Figure  5.  Views  of  computational  mesh:  (a)  tail  region  of  grid 
fin  case  and  (b)  close-up  of  grid  fin  base. 

Mesh  adaption  was  used  to  determine  mesh  independence.  The  meshes  for  two  planar  fin  cases 
(a  =  4°  and  10°;  both  at  8  =  10°)  were  adapted  on  the  static  pressure  gradient  using  the  adaption 
tools  in  the  FLUENT  software.  After  adaption,  the  solutions  were  again  run  until  convergence 
was  achieved.  The  maximum  change  in  the  aerodynamic  coefficients  was  0.9%  for  the  a  =  4° 
case,  with  four  of  the  six  coefficients  changing  less  than  0.5%.  The  coefficients  changed  less 
than  0.2%  after  mesh  adaption  in  the  a  =  10°  case.  These  results  show  that  the  original  mesh 
used  for  the  planar  fin  case  had  high  enough  resolution  for  a  mesh  independent  solution.  For  the 
grid  fin  case,  the  mesh  on  the  first  13  cal.  is  the  same  as  the  planar  fin  case,  so  it  is  also  adequate. 
In  the  tail  region,  the  small  cell  sizes  required  for  the  grid  fin  surfaces  resulted  in  a  volume  mesh 
that  should  have  had  more  than  adequate  resolution,  so  no  mesh  adaption  was  performed. 
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3.  Results  and  Discussion 


3.1  Aerodynamic  Coefficients  -  CFD  Validation 

The  viscous  and  pressure  forces  from  the  computed  flowfields  were  integrated  along  the  missile 
body,  canard,  and  fin  surfaces  to  calculate  the  aerodynamic  coefficients.  The  axial  force,  normal 
force,  side  force,  rolling  moment,  pitching  moment,  and  yawing  moment  are  presented  in 
missile-based  coordinates.  This  is  a  right-handed  system  with  the  x-axis  coinciding  with  the 
missile  axis  and  oriented  to  the  rear,  the y-axis  oriented  to  the  missile’s  starboard  side,  and  the  z- 
axis  oriented  upward.  The  forces  are  positive  when  coinciding  with  the  positive  coordinate  axes. 
By  convention,  the  moments  are  defined  as  follows:  the  rolling  moment  is  positive  when  the  roll 
is  clockwise  when  looking  forward  from  the  aft  end  of  the  missile,  the  yawing  moment  is 
positive  when  the  nose  is  moving  right,  and  the  pitching  moment  is  positive  when  the  nose  is 
moving  upward.  The  reference  area  is  the  cross-sectional  area  of  the  missile  base,  and  the 
reference  length  is  the  diameter  of  the  missile.  The  calculated  coefficients  are  compared  to 
wind-tunnel  measurements  performed  at  DREV  [5],  Since  the  base  flow  was  not  calculated, 
only  the  “forebody”  axial  force  is  compared  to  the  data.  The  forebody  axial  force  from  the 
experimental  data  was  obtained  by  correcting  the  measured  total  axial  force  for  the  base  drag 
component  using  measured  pressures  at  the  base  of  the  wind  tunnel  model.  Not  including  the 
base  flow  in  the  simulations  can  result  in  error  of  a  few  percent  in  the  computed  forebody  axial 
force,  while  using  localized  base  pressure  measurements  to  obtain  the  experimental  forebody 
axial  force  can  also  lead  to  errors  of  a  few  percent  [10]. 

All  the  computed  aerodynamic  coefficients  showed  very  good  agreement  with  experimental 
values  from  the  DREV  wind  tunnel  (Appendices  A  and  B).  The  computed  and  measured  normal 
force  and  pitching  moment  coefficients  are  compared  in  Figure  6  and  Figure  7,  respectively,  for 
theM=  1.5  and  M=  3.0  planar  fin  cases.  Also  shown  in  Figures  6  and  7  are  the  predictions  from 
the  AP98  aeroprediction  code  [13].  AP98  is  an  engineering  code  based  on  empirical  and 
theoretical  methods.  The  comparison  of  the  CFD  and  AP98  predictions  is  excellent.  The 
difference  between  the  predicted  and  measured  values  increases  with  a.  In  general,  the  slope  of 
the  linear  portion  of  the  curve  is  steeper  in  the  predicted  case.  However,  the  predicted  trend  is 
the  same  as  the  measured  trend,  including  the  increased  curvature  at  the  higher  Mach  number. 
Similar  results  were  obtained  for  the  grid  fin  cases,  which  are  presented  in  Appendix  B. 

The  computed  and  measured  axial  force  coefficients  are  compared  in  Figure  8  and  Figure  9  for 
theM=  1.5  planar  and  grid  fin  cases,  respectively.  Again  the  CFD  and  AP98  predictions 
compare  very  well.  The  difference  between  the  predicted  and  measured  Cx  values  ranged  from 
about  7.5%-15%.  The  forebody  Cx  for  the  grid  fin  case  is  about  1.7  times  that  observed  in  the 
planar  fin  case  at  M=  1 .5.  Note  that  the  drag  of  the  grid  fins  on  this  model  is  larger  than  that 
which  could  be  achieved  in  an  “optimum”  design.  First,  the  grid  fin  frame  is  not  chamfered, 
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(a) 


(b) 


Figure  6.  Normal  force  coefficient,  planar  fin  case,  (a)  M  =  1.5  and  (b)  M  =  3.0. 

which  has  been  shown  to  dramatically  reduce  the  fin  drag  [8].  Second,  because  of  the  small  size 
of  the  wind  tunnel  model,  the  web  thickness  could  not  be  scaled  down  to  the  proper  design 
thickness  due  to  machining  limitations,  and  it  is  about  1.5  times  larger  than  optimum.  The  fin 
drag  was  previously  shown  to  be  directly  related  to  the  thickness  of  the  webs  and  the  number  of 
webs  present  [14].  Earlier  studies  [9, 10],  with  an  optimum  web  thickness  and  chamfered  grid 
fin  frame,  showed  an  increase  in  forebody  Cx  for  a  missile  with  grid  fins  of  only  1 .4  times  that  of 
a  missile  with  planar  fins. 
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(a) 


(b) 

Figure  7.  Pitching  moment  coefficient  about  nose,  planar  fin  case,  (a)  M=  1.5  and  (b) 
M=  3.0. 
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Figure  9.  Axial  force  coefficient,  grid  fin  case,  M =  1.5. 


An  important  objective  of  this  study  was  to  accurately  predict  the  adverse  induced  rolling 
moment  and  side  force  observed  in  the  wind  tunnel  experiments.  Validation  of  these 
components  gives  confidence  to  the  flow  visualizations  obtained  from  the  CFD.  As  shown  in 
Figure  10,  very  good  agreement  between  the  computed  and  experimental  values  of  Q  was 
obtained.  Figure  1 0  shows  the  comparison  of  the  computed  and  measured  Q  for  the  planar  and 
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(b) 

Figure  10.  Rolling  moment  coefficient  for  (a)  planar  fin  and  (b)  grid  fin  cases. 
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grid  fin  cases  at  both  Mach  numbers  and  5  =  10°.  The  experimental  Q  had  a  large  amount  of 
bias.  This  bias  was  removed  for  the  8=10°  case  shown  in  Figure  10,  using  a  method  that  DREV 
uses,  which  was  as  follows.  The  Q  for  the  8  =  0°  case,  which  should  be  zero  at  all  a,  is 
subtracted  from  the  Q  for  the  8  =  10°  case.  This  assumes  that  the  bias  is  constant  between  test 
runs.  A  small  difference  between  the  computed  and  experimental  Q  remains;  however,  note  that 
the  absolute  magnitude  of  the  Q  is  small  (<0.1),  so  the  agreement  is  still  quite  good.  The  CFD 
accurately  captured  the  adverse  roll  phenomenon.  The  agreement  in  the  shape  of  the  curve, 
including  reversal  of  the  roll  direction  in  the  M=  1 .5  planar  fin  case  (Figure  10a),  is  excellent. 
Also  note  that  the  direction  of  roll  in  the  M—  1.5  planar  fin  case  is  opposite  that  which  should  be 
induced  by  the  deflected  canards  until  a  is  above  8°.  This  effect,  which  appears  to  be  due  to  the 
canard  trailing  vortices  interacting  with  the  planar  fins,  will  be  discussed  further  when  the 
flowfield  is  presented.  For  the  M  =  1.5  grid  fin  case  (Figure  10b),  Q  is  positive  at  a  =  0°, 
decreasing  to  near  zero  between  4°  <  a  <  7°  (some  loss  of  roll  control).  At  M =  3.0,  Q  is  always 
positive  and  is  similar  for  both  fin  types. 

The  CFD  also  accurately  predicted  the  nonlinear  induced  side  force  characteristics  observed  in 
the  canard  deflected  cases.  Figure  1 1  shows  the  comparison  of  the  computed  and  measured  Cy 
for  the  planar  and  grid  fin  cases  at  both  Mach  numbers  and  8  =  10°.  The  flow  is  symmetric 
about  the  vertical  plane  for  the  8  =  0°  cases,  so  Cy  =  0  and  is  not  shown.  There  was  little 
difference  in  Cy  due  to  fin  type  at  either  Mach  number.  At  M=  3.0,  the  maximum  Cy  was  not  too 
large.  However,  atM=  1.5,  Cy  was  about  5  times  larger  than  at  M- 3.0.  The  calculations 
accurately  predicted  the  change  in  slope  of  the  Cy  vs.  a  curve  at  M  =  3.0  for  both  fin  types.  The 
predictions  do  not  predict  maximum  Cy  as  well  at  M-  1 .5.  The  simulations  end  near  the 
maximum  value  but  it  appears  that  the  beginning  of  the  change  in  slope  is  correctly  predicted  in 
the  grid  fin  case  (Figure  1  lb). 

The  comparisons  of  the  yawing  moment  were  similar  to  the  side  force  and  are  shown  in 
Appendices  A  and  B,  where  all  coefficients  for  both  fin  types  are  presented. 

These  results  show  that  the  CFD  has  accurately  predicted  the  aerodynamic  coefficients  for  a 
generic  canard-controlled  missile  with  both  conventional  planar  fins  and  grid  fins  in  supersonic 
flow.  With  the  CFD  results  validated,  the  computed  flowfield  can  now  be  used  to  investigate  the 
physics  responsible  for  roll-reversal  or  loss  of  roll  control. 

3.2  Flowfield  Visualizations 

Visualizations  of  the  flowfield  showed  that  the  canard  deflection  had  a  strong  effect  on  the  forces 
on  the  missile  body.  Figure  12  shows  the  Cp  distribution  on  the  missile  surfaces  for  the  planar 
fin  cases  with  8  =  10°  at  a  =  4°  and  10°.  The  8  =  0°  cases  are  not  presented,  but  they  showed  a 
periodic  Cp  distribution  on  the  missile  with  no  contribution  to  Cy,  Ci,  or  C„.  For  the  8  =  10° 
cases,  the  canard  wake  and  trailing  vortices  interacted  with  the  missile  flowfield  to  change  the 
pressure  distribution  both  along  the  missile  body  and  on  the  tail  fins.  The  effect  is  most 
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Figure  11.  Side  force  coefficient  for  (a)  planar  fin  and  (b)  grid  fin  cases. 
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Figure  12.  Port  (left  image)  and  starboard  (right  image)  views  of  Cp  contours  at 

M  -  1.5  for  planar  fin  case  with  8  =  10°  and  (a,  b)  a  =  4°  and  (c,  d)  a  =  10°. 

pronounced  at  M  =  1 .5,  where  a  large  low-pressure  region  is  observed  on  the  starboard  side  of 
the  missile  (Figure  12d).  Although  the  effect  was  still  present  at  M  =  3.0,  the  effect  was  not  as 
great,  which  correlates  with  the  measured  Cy. 

The  large  side  force  at  M  =  1 .5  and  high  a  appears  to  be  caused  by  the  interaction  of  the  wake 
flow  from  the  deflected  canards  with  the  crossflow  separation  on  the  leeward  side  of  the  missile. 
This  is  best  observed  by  viewing  visualizations  of  the  pressure  distribution  on  an  axial  cut-plane 
being  swept  along  the  missile  axis.  The  plots  in  Figure  13  are  images  from  such  a  visualization 
and  also  illustrate  the  effect.  As  the  angle  of  attack  of  the  missile  increases,  the  effective  angle 
of  attack  of  the  starboard  canard  decreases  while  that  of  the  port  canard  increases.  The  wake 
from  the  port  canard  appears  to  be  primarily  responsible  for  inducing  a  clockwise  (viewed  from 
rear  of  missile)  flow  around  the  missile  behind  the  canards  that  both  shifts  the  azimuthal  location 
and  increases  the  size  of  the  low-pressure  region  that  is  due  to  the  crossflow  around  the  missile. 
The  leeward  canard  also  appears  to  contribute  to  this  effect.  The  low-pressure  wake  behind  both 
the  port  and  leeward  canards  increase  as  a  increases  (Figure  13a-c). 
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Figure  13.  Cp  contours  on  missile  surface  and  axial  plane  behind  canards 
for  8  =  10°  and  a  =  0°,  4°,  and  10°  at  (a-c)  M=  1.5  and  (d-f) 
Af  =  3.0. 


The  canard  deflection  should  give  the  missile  a  clockwise  (when  viewed  from  rear),  or  positive, 
rolling  moment.  At  M  =  1 .5,  the  Q  is  actually  negative  until  a  >  8°.  This  adverse  roll  can  be 
explained  by  examining  the  pressure  distribution  and  forces  on  the  tail  fins.  A  high-pressure 
region  on  the  starboard  side  of  the  leeward  fin  at  a  =  4°  can  be  seen  in  Figure  14a  (and 
Figure  12b).  Although  not  shown  in  the  figures,  a  similar  high-pressure  region  was  also  on  the 
windward  side  of  the  starboard  fin.  These  unbalanced  pressures  on  the  tail  fins  are  enough  to 
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Figure  14.  Cp  contours  on  planar  fins  at  M  =  1.5  and  at  (a) 
a  =  4°  and  (b)  a  =  10°. 


reverse  the  rolling  moment.  At  a  =  10°,  the  effect  is  much  lower,  which  appears  to  be  primarily 
due  to  the  vortices  trailing  the  port,  starboard,  and  leeward  canards  missing  the  tail  fins  entirely. 
In  fact,  there  are  unbalanced  forces  (Figure  14b  and  Figure  12c)  that  act  in  addition  to  the  canard 
forces  to  induce  added  positive  roll. 

Similar  flowfield  results  were  observed  for  the  grid  fin  case  (Figure  15).  Although  there  was  still 
a  nonuniform  pressure  distribution  on  the  fins,  the  characteristic  geometry  of  the  grid  fins 
alleviated  the  unwanted  rolling  moment  effect.  At  a  =  0,  the  unwanted  rolling  moment  at 
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(b) 

Figure  15.  Cp  contours  on  grid  fins  at  M  =1.5  and  at  (a) 
a  =  4°  and  (b)  a  =  10°. 


M  =  1.5  observed  in  the  planar  fin  case  is  eliminated  (Figure  10).  There  is  still  a  range  between 
4°  <  a  <  7°  that  the  effectiveness  of  the  canards  is  severely  reduced,  but  there  is  still  an 
improvement  over  the  planar  tail  fins.  These  results  are  shown  quantitatively  in  section  3.3, 
where  the  forces  on  the  control  surfaces  are  presented. 

Scientific  visualization  movies  indicated  that  the  adverse  forces  on  the  missile  and  fins  were  due 
to  a  combination  of  the  canard-trailing  vortices  interacting  with  both  the  tail  fins  and  the 
crossflow  separation  region  on  the  leeward  side  of  the  missile.  The  lack  of  adverse  rolling 
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moment  effects  at  M  =  3.0  appears  to  be  due  directly  to  the  canard-trailing  vortices  having  a 
lower  strength  at  the  higher  Mach  number.  Plots  of  vorticity  magnitude  at  several  axial  planes 
are  shown  in  Figure  16  for  the  M=  1.5  cases  and  Figure  17  for  the  M  =  3.0  cases,  both  with 
6  =  10°.  At  M  =  1.5  (Figure  16)  the  canard-trailing  vortices  have  retained  significant  strength 
when  they  reach  the  tail  fins.  However,  at  M  =  3.0  (Figure  17)  the  canard-trailing  vortices  have 
lost  much  of  their  strength  by  the  time  they  reach  halfway  down  the  missile  length  (x/D  =  8). 


Figure  16.  Cp contours  (range:  -0.1-0.3)  on  missile  surfaces  and  vorticity  magnitude  contours  (range: 

0-50000)  on  axial  cross-planes  at  2,  4, 6,  8, 10, 12, 14,  and  16  cal.  for  planar  fin  case  with 
8=  10°  at  M  =  1.5  and  (a)  a  =  0°,  (b)  a  =  4°,  and  (c)  a=  10°. 


Other  than  the  lower  vorticity  magnitude  at  the  higher  Mach  number,  the  general  characteristics 
of  the  trailing  vortices  are  similar  at  both  Mach  numbers.  At  a  =  0°,  the  trailing  vortices  shift 
counterclockwise  (when  viewed  from  rear)  and  pass  between  the  fins,  but  closer  to  the  fin 
directly  behind  the  respective  canard.  There  also  appears  to  be  a  small  interaction  with  the 
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Figure  17.  Cp  contours  (range:  -0. 1-0.3)  on  missile  surfaces  and  vorticity  magnitude  contours  (range: 

0-50000)  on  axial  cross-planes  at  2, 4,  6,  8,  10,  12,  14,  and  16  cal.  for  planar  fin  case  with 
6  =  10°  at  M  =  3.0  and  (a)  a  =  0°,  (b)  a  =  4°,  and  (c)  a  =  10°. 

missile  boundary  layer.  At  a  =  4°,  the  windward  canard-trailing  vortex  merges  with  the 
boundary  layer  vorticity  at  about  x/D  =  8.  The  leeward  canard-trailing  vortex  passes  near  the  tip 
of  the  leeward  fin.  The  trailing  vortices  from  the  port  and  starboard  canards  move  to  the  leeward 
side  of  the  missile  and  begin  to  merge  with  the  leeward  canard-trailing  vortex  just  behind  the 
leeward  fin.  At  a  =  10°,  the  windward  canard-trailing  vortex  merges  with  the  boundary  layer 
vorticity  at  about  x/D  =  4.  The  vortices  from  the  other  three  canards  merge  on  the  leeward  side 
of  the  missile  at  about  x/D  =  10  and  pass  over  the  tail  fins.  At  a  =  10°,  the  crossflow  separation 
vortices  primarily  interact  with  the  leeward  tail  fin.  From  these  visualizations,  it  appears  that  the 
interaction  of  the  trailing  vortices  with  the  tail  fins  is  primarily  responsible  for  the  adverse  forces 
on  the  tail  fins  that  affect  the  rolling  moment. 

For  the  8  =  0°  cases  (Figure  18)  at  M  =  1.5,  the  strength  of  the  canard-trailing  vortices  is  less 
than  that  observed  for  the  8=10°  cases.  The  strength  of  the  port  and  starboard  canard-trailing 
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Figure  18.  Cp  contours  (range:  -0. 1-0.3)  on  missile  surfaces  and  vorticity  magnitude 
contours  (range:  0-50000)  on  axial  cross-planes  at  2, 4,  6,  8,  10, 12, 14,  and 
16  cal.  for  planar  fin  case  with  8  =  0°  at  M  =  1 .5  and  (a)  a  =  0°  and  (b)  a  =  4°. 

vortices  increases  as  a  increases;  however,  the  flow  is  symmetric  near  the  tail  fins  and  no 
adverse  rolling  moment  is  produced. 

The  effect  of  the  grid  fins  on  the  canard-trailing  vortices  is  shown  in  Figure  19,  which  compares 
the  vorticity  magnitude  on  several  cross-planes  near  the  tail  of  the  missile  for  the  planar  and  grid 
fin  cases  at  M  =  1.5,  a  =  4°,  and  8  =  10°.  In  the  planar  fin  case  (Figure  19a)  the  trailing  vortices 
from  the  port,  starboard,  and  leeward  fins  are  nearly  unchanged  after  they  pass  by  the  leeward 
tail  fin.  However,  in  the  grid  fin  case  (Figure  19b),  although  the  trailing  vortex  from  the  port  fin 
passes  between  the  port  and  leeward  fin  relatively  unchanged,  the  other  two  vortices  are 
dramatically  changed  from  the  interaction  with  the  leeward  grid  fin.  The  force  on  the  fin  from 
this  interaction  does  not  favor  any  particular  direction,  so  a  large  moment  is  not  produced. 
However,  in  the  planar  fin  case,  the  localized  high  pressure  on  the  leeward  fin  does  cause  the 
adverse  rolling  moment. 

3.3  Control  Surface  Forces 

The  forces  on  the  canards  and  fins  are  summarized  in  tabular  and  graphical  form  in 
Appendices  C-F.  Some  results  are  presented  in  this  section  to  quantify  the  trends  observed  in 
the  flow  visualizations. 

The  forces  on  the  leeward  and  windward  fins  resulting  in  the  adverse  rolling  moment  are  shown 
in  Figure  20,  which  shows  the  side  force  coefficient  on  each  fin  for  the  planar  and  grid  fin  cases 
at  8  =  10°.  In  the  planar  fin  case,  the  windward  fin  (fin2)  is  at  a  nearly  constant,  small  positive 
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Figure  19.  Cp  contours  (range:  -0.1-0.3)  on  missile  surfaces  and  vorticity 

magnitude  contours  (range:  0-50000)  on  axial  cross-planes  at  (a)  12, 
14,  and  16  cal.  for  the  planar  fin  case  and  at  (b)  12, 14,  and  15  cal.  for 
the  grid  fin  case  with  8  =  10°  at  M  =  1 .5  and  a  =  4°. 


(b) 

Figure  20.  Side  force  coefficient  on  tail  fins  for  (a)  planar  fin  and  (b)  grid  fin  case  for 
8  =  10°  and  M=  1.5. 
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value.  The  leeward  fin  (fin  4)  is  negative  until  a  >  8°.  This  force  imbalance  clearly  produces  the 
adverse  roll  since  the  fin  force  coefficients  are  of  the  same  order  as  those  on  the  canards  and  the 
larger  moment  arm  of  the  fins  therefore  produces  moments  larger  than  those  produced  by  the 
canards.  In  the  grid  fin  case  (Figure  20b),  there  are  still  side  forces  on  the  windward  and  leeward 
fins  that  tend  to  induce  a  negative  roll,  but  this  time  they  are  lower  and  do  not  reduce  the  roll 
effectiveness  of  the  canards  as  much.  Also,  observe  that  the  horizontal  grid  fins  produce  a  side 
force,  unlike  the  planar  fins.  This  does  not  affect  the  roll  moment  and  only  adds  a  very  small 
contribution  to  missile  side  force. 

Similarly,  but  to  a  smaller  extent,  the  normal  force  on  the  starboard  and  port  planar  fins 
contributes  to  the  adverse  rolling  moment.  This  is  illustrated  in  Figure  21,  which  shows  that  the 
difference  in  normal  force  between  the  port  and  starboard  grid  fins  is  less  than  that  for  the  planar 
fins. 

One  advantage  noted  for  grid  fins  over  planar  fins  is  that  they  maintain  their  effectiveness  at  high 
Mach  number.  This  can  be  observed  by  comparing  Figure  21  and  Figure  22,  which  show  the  fin 
normal  force  for  both  fin  types  and  both  Mach  numbers.  For  this  particular  grid  fin  design,  at 
a  =  10°  and  M=  1.5,  the  total  fin  normal  force  in  the  grid  fin  case  is  about  65%  that  in  the  planar 
fin  case.  At  M-  3.0,  the  total  fin  normal  force  in  the  planar  fin  case  has  decreased  over  50%. 
However,  the  total  fin  normal  force  in  the  grid  fin  case  has  stayed  nearly  the  same  as  that  at 
M=  1.5.  The  small  reduction  is  mainly  due  to  the  normal  force  on  the  leeward  fin  going 
negative  at  higher  a. 


4.  Summary  and  Conclusions 


Viscous  CFD  calculations  were  used  to  predict  the  aerodynamic  coefficients  and  flowfield 
around  a  generic  canard-controlled  missile  configuration  in  supersonic  flow.  Validation  of  the 
computed  results  was  demonstrated  by  the  very  good  agreement  between  the  computed 
aerodynamic  coefficients  and  those  obtained  from  wind  tunnel  measurements. 

Visualizations  of  the  flowfield  at  low  supersonic  speed  showed  that  the  canard  downwash  and 
trailing  vortices  produced  a  low-pressure  region  on  the  starboard  side  of  the  missile  that,  in  turn, 
produced  an  adverse  induced  side  force.  The  side  force  was  much  lower  at  the  higher  supersonic 
speed. 

The  visualizations  also  showed  that  the  canard-trailing  vortices  interact  with  the  tail  fins  until  a 
is  high  enough  so  that  the  vortices  miss  the  leeward  fin.  The  pressure  differential  on  the  leeward 
fin  produced  by  the  interaction  with  the  canard-trailing  vortices  is  primarily  responsible  for  the 
adverse  induced  roll  effects  observed  when  planar  fins  are  used. 
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(b) 

Figure  21.  Normal  force  coefficient  on  tail  fins  for  (a)  planar  fin  and  (b)  grid  fin  case 
5  =  10°  and  M=  1.5. 
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Grid  tail  fins  improved  the  roll-control  effectiveness  of  the  canards  at  low  supersonic  speed. 
There  was  still  an  uneven  pressure  distribution  on  the  grid  tail  fins,  produced  by  the  interaction 
with  the  canard-trailing  vortices.  However,  the  design  of  the  grid  fin  distributes  the  pressure 
differently  so  that  the  side  forces  are  lower  and  do  not  produce  as  large  a  roll  moment  as  the 
planar  fins  do. 

No  adverse  rolling  moment  was  observed  with  no  canard  deflection  or  at  the  higher  supersonic 
speed  for  either  tail  fin  type.  Flow  visualizations  showed  that  the  intensity  of  the  canard-trailing 
vortices  was  much  lower  in  these  two  cases.  Methods  to  reduce  the  intensity  of  the  canard 
trailing  vortices,  or  direct  them  away  from  the  tail  fins,  may  be  two  possible  ways  to  alleviate  the 
roll-control  problems  in  munitions  using  canards  for  control. 

The  validated  CFD  results  confirm  that  the  use  of  grid  tail  fins  dramatically  improves  the  roll 
effectiveness  of  the  canards  at  low  supersonic  speed.  Visualizations  of  the  flowfield  from  the 
simulations  performed  in  this  study  help  in  the  understanding  of  the  flow  physics  and  can  lead  to 
improved  tail  fin  and  canard  designs. 
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Appendix  A.  Aerodynamic  Coefficients  for  Planar  Fin  Case 
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Figure  A-L  Normal  force  for  planar  fin  case  at  Mach  1.5. 


Figure  A-2.  Normal  force  for  planar  fin  case  at  Mach  3.0. 
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Figure  A-6.  Axial  force  for  planar  fin  case  at  Mach  3.0. 
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Figure  A-7.  Pitching  moment  about  nose  for  planar  fin  case  at  Mach  1.5. 


Figure  A-8.  Pitching  moment  about  nose  for  planar  fin  case  at  Mach  3.0. 
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Figure  A-9.  Pitching  moment  about  moment  reference  point  (MRP)  for  planar  fin  case  at 
Mach  1.5. 


Figure  A- 10.  Pitching  moment  about  MRP  for  planar  fin  case  at  Mach  3.0. 
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Figure  A-13.  Yawing  moment  for  planar  fin  case  at  Mach  1 .5. 


Figure  A-14.  Yawing  moment  for  planar  fin  case  at  Mach  3.0. 
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Figure  A-16.  Center  of  pressure  location  from  nose  for  planar  fin  case  at  Mach  3.0. 
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Appendix  B.  Aerodynamic  Coefficients  for  Grid  Fin  Case 
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Figure  B-l .  Normal  force  for  grid  fin  case  at  Mach  1 .5. 


Figure  B-2.  Normal  force  for  grid  fin  case  at  Mach  3.0. 
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Figure  B-3.  Side  force  for  grid  fin  case  at  Mach  1.5. 


Figure  B-4.  Side  force  for  grid  fin  case  at  Mach  3.0. 
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Figure  B-6.  Axial  force  for  grid  fin  case  at  Mach  3.0. 
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Figure  B-8.  Pitching  moment  about  nose  for  grid  fin  case  at  Mach  3.0. 
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Figure  B-9.  Pitching  moment  about  moment  reference  point  (MRP)  for  grid  fin  case  at 
Mach  1.5. 


Figure  B-10.  Pitching  moment  about  MRP  for  grid  fin  case  at  Mach  3.0. 
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DREV  EXP  (M=3.0,  Canard=0°) 
-B*  *  CFD  (M=3  0,  Canard=0°) 

DREV  EXP  (M=3  0,  Canard=10° 
-A-  CFD  (M=3  0,  Canard=10°) 


Figure  B-14.  Yawing  moment  for  grid  fin  case  at  Mach  3.0. 
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Figure  B-15.  Center  of  pressure  location  from  nose  for  grid  fin  case  at  Mach  1.5. 
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Appendix  C.  Force  Coefficients  on  Canards 
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Figure  C-l.  Canard  normal  force,  8  =  0°,  Mach  1.5. 


Figure  C-2.  Canard  normal  force,  8  =  0°,  Mach  3.0. 


45 


46 


47 


0.08 

0  06 
004 
0  02 


Cy  0  00 


-0  02 


-0  04 
-0  06 


-0  08 


-canardl  (S) 
-canard2  (W) 
-canard3  (P) 
-canard4  (L) 

■  Total  Canard 


6 

a 


Figure  C-7.  Canard  side  force,  5  =  10°,  Mach  1.5. 


Figure  C-8.  Canard  side  force,  5  =  10°,  Mach  3.0. 
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Figure  C-9.  Canard  axial  force,  8  =  0°,  Mach  1.5. 


Figure  C-10.  Canard  axial  force,  8  =  0°,  Mach  3.0. 
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Figure  C-12.  Canard  axial  force,  8  =  10°,  Mach  3.0. 
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Table  C-l.  Aerodynamic  coefficients  on  canards,  planar  fin  case,  6  =  0°,  Mach  1.5. 


a  0 

2 

4 

6 

8 

10 

(Axial  Force 

Canard  1  (S) 

■nmkHfc 

0.0062 

0.0062 

0.0062 

0.0061 

0.0059 

Canard2  (W) 

PPi 

0.0062 

0.0064 

0.0065 

0.0066 

0.0067 

Canard3  (P) 

0.0061 

0.0062 

0.0062 

0.0062 

0.0061 

0.0059 

Canard4  (L) 

0.0061 

0.0060 

0.0059 

0.0058 

0.0056 

0.0055 

Total  canard 

0.0245 

0.0246 

0.0247 

0.0246 

0.0244 

0.0241 

Side  Force 

Canard  1  (S) 

-0.0007 

-0.0007 

-0.0006 

-0.0004 

-0.0002 

-0.0001 

Canard2  (W) 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Canard3  (P) 

0.0007 

0.0007 

0.0006 

0.0004 

0.0002 

0.0001 

Canard4  (L) 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Total  canard 

0.0000 

0.0000 

■ittlllllllM 

mm 

Normal  Force 

Canardl  (S) 

0.0000 

0.0173 

0.0341 

0.0501 

0.0650 

0.0790 

Canard2  (W) 

0.0007 

0.0009 

0.0011 

0.0014 

0.0016 

0.0019 

Canard3  (P) 

0.0000 

0.0173 

0.0341 

0.0501 

0.0650 

0.0789 

Canard4  (L) 

-0.0007 

-0.0004 

-0.0002 

-0.0001 

EEMH 

Total  canard 

mmm 

msEsm 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 


Table  C-2.  Aerodynamic  coefficients  on  canards,  planar  fin  case,  8  =  10°,  Mach  1.5. 


a  0 

2 

4 

6  8 

10 

(Axial  Force 

Canardl  (S) 

0.0159 

0.0132 

0.0101 

0.0067 

0.0032 

0.0002 

Canard2(W) 

0.0159 

0.0162 

0.0164 

0.0166 

0.0168 

0.0169 

Canard3  (P) 

0.0159 

0.0183 

0.0204 

0.0222 

0.0239 

0.0254 

Canard4  (L) 

0.0159 

0.0157 

0.0152 

0.0151 

0.0149 

Total  canard 

WiiiT&y 

■HI 

KSiaZI 

Side  Force 

Canardl  (S) 

-0.0002 

-0.0004 

-0.0006 

-0.0007 

-0.0009 

-0.0009 

Canard2  (W) 

-0.0548 

-0.0556 

-0.0564 

-0.0572 

-0.0578 

-0.0584 

0.0001 

-0.0001 

-0.0002 

-0.0002 

-0.0003 

Canard4  (L) 

0.0540 

0.0534 

0.0530 

0.0527 

0.0525 

Total  canard 

0.0000 

-0.0019 

-0.0036 

-0.0051 

-0.0062 

-0.0071 

Normal  Force 

Canardl  (S) 

-0.0548 

-0.0393 

-0.0224 

-0.0041 

0.0155 

0.0350 

Canard2  (W) 

0.0002 

0.0004 

0.0005 

0.0007 

0.0009 

0.0010 

Canard3  (P) 

0.0548 

0.0689 

0.0818 

0.0937 

0.1049 

0.1155 

Canard4  (L) 

-0.0002 

-0.0001 

0.0000 

0.0001 

0.0002 

0.0003 

Total  canard 

0.0000 

0.0298 

0.0599 

0.0904 

0.1215 

0.1519 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward 
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Table  C-3.  Aerodynamic  coefficients  on  canards,  planar  fin  case,  5  =  0°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Canard  1  (S) 

0.0058 

0.0058 

0.0058 

0.0059 

IBSEESI 

0.0061 

Canard2  (W) 

0.0058 

0.0061 

0.0064 

0.0066 

Canard3  (P) 

0.0058 

0.0058 

0.0058 

0.0059 

0.0060 

ifiireMi 

Canard4  (L) 

0.0051 

0.0048 

0.0045 

life'll 

Total  canard 

HI 

0.0232 

■m 

Side  Force 

Canardl  (S) 

-0.0010 

-0.0010 

-0.0010 

-0.0011 

-0.0011 

-0.0012 

Canard2  (W) 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Canard3  (P) 

0.0010 

0.0010 

0.0010 

0.0011 

0.0011 

0.0012 

Canard4  (L) 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Total  canard 

0.0000 

0.0000 

HlklllllllM 

rflRlIHlIlM 

Normal  Force 

Canardl  (S) 

0.0000 

0.0074 

0.0148 

0.0222 

0.0299 

0.0378 

Canard2  (W) 

0.0010 

0.0012 

0.0013 

0.0017 

0.0018 

Canard3  (P) 

0.0000 

0.0074 

0.0148 

0.0299 

0.0378 

Canard4  (L) 

-0.0010 

-0.0008 

-0.0006 

-0.0001 

Total  canard 

HI 

0.0303 

mmm 

0.0773 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 


Table  C-4.  Aerodynamic  coefficients  on  canards,  planar  fin  case,  5=10°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Canardl  (S) 

0.0117 

0.0104 

0.0091 

0.0066 

0.0055 

Canard2  (W) 

0.0117 

0.0126 

0.0133 

0.0146 

0.0152 

Canard3  (P) 

0.0118 

0.0132 

0.0148 

0.0164 

0.0182 

0.0200 

Canard4  (L) 

0.0117 

0.0109 

0.0102 

0.0095 

0.0089 

0.0083 

Total  canard 

filHM 

0.0471 

0.0473 

0.0477 

0.0483 

Side  Force 

Canardl  (S) 

besom 

-0.0011 

f  ifflfM 

-0.0010 

-0.0010 

Canard2  (W) 

PIXIXW 

-0.0347 

-0.0420 

-0.0443 

Canard3  (P) 

0.0011 

final 

0.0013 

0.0014 

0.0016 

0.0017 

Canard4  (L) 

0.0321 

0.0274 

0.0254 

0.0236 

0.0220 

Total  canard 

0.0000 

-0.0050 

-0.0096 

-0.0139 

-0.0178 

-0.0215 

Normal  Force 


Canardl  (S) 

-0.0321 

-0.0244 

-0.0171 

-0.0100 

0.0039 

Canard2  (W) 

0.0011 

0.0013 

0.0015 

0.0017 

0.0021 

Canardl  (P) 

0.0321 

0.0400 

0.0484 

f  iEftM 

0.0750 

Canard4  (L) 

-0.0011 

-0.0009 

-0.0007 

B3BS3I 

PE2SS3I 

-0.0002 

Total  canard 

0.0000 

0.0160 

0.0320 

0.0481 

0.0644 

0.0807 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 
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Appendix  D.  Force  Coefficients  on  Planar  Fins 


Figure  D-l.  Planar  fin  normal  force,  5  =  0°,  Mach  1.5. 


Figure  D-2.  Planar  fin  normal  force,  8  =  0°,  Mach  3.0. 
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Figure  D-4.  Planar  fin  normal  force,  8  =  10°,  Mach  3.0 
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Figure  D-8  Planar  fin  side  force,  8  =  10°,  Mach  3.0 
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Figure  D-10.  Planar  fin  axial  force,  8  =  0°,  Mach  3.0. 
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Figure  D-l  1 .  Planar  fin  axial  force,  5  =  10°,  Mach  1 .5. 
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Figure  D-12.  Planar  fin  axial  force,  5  =  10°,  Mach  3.0. 
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Table  D-l.  Aerodynamic  coefficients  on  planar  fins,  5  =  0°,  Mach  1.5. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Finl  (S) 

0.0207 

0.0208 

IBifMB 

0.0227 

0.0229 

Fin2(W) 

0.0207 

0.0211 

0.0213 

0.0215 

0.0215 

0.0214 

Fin3  (P) 

0.0207 

0.0208 

0.0215 

0  0221 

0.0227 

0.0229 

Fin4  (L) 

0.0207 

0.0199 

0.0202 

0.0206 

0.0208 

0.0208 

Total  fin 

0.0826 

0.0828 

0.0844 

0.0864 

0.0878 

0.0880 

Side  Force 

Finl  (S) 

0.0000 

0.0005 

mm i 

0.0032 

0.0057 

Fin2  (W) 

0.0000 

0.0000 

mmm 

Fin3  (P) 

0.0000 

-0.0005 

-0  0017 

-0.0046 

BBI 

Fin4  (L) 

0.0000 

0.0000 

0  0000 

-0  0001 

-0.0001 

-0.0002 

Total  fin 

0.0000 

0.0000 

0.0000 

-0.0001 

-0.0001 

-0.0002 

Normal  Force 

Finl  (S) 

0.0000 

0.0832 

01712 

Fin2  (W) 

0  0007 

Fin3  (P) 

0.0832 

PH 

Fin4  (L) 

0  0002 

0  0001 

-0.0001 

Total  fin 

0.0000 

0.1670 

0.3433 

0.5228 

0.6972 

0.8575 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 


Table  D-2.  Aerodynamic  coefficients  on  planar  fins,  5  =  10°,  Mach  1.5. 


a  0 

2 

4 

6  8 

10 

Axial  Force 

Finl  (S) 

0.0212 

0.0210 

0.0225 

Fin2  (W) 

0  0211 

0.0212 

IMI 

0.0215 

mmm 

Mmm 

Fin3  (P) 

0.0211 

0.0205 

0.0210 

0.0218 

0  0224 

0.0227 

Fin4  (L) 

B031M 

BCTilWSi 

0.0188 

0.0189 

Total  fin 

WUZZM 

0.0833 

'mi 

Side  Force 

Finl  (S) 

-0.0001 

WM2MM 

Fin2  (W) 

0.0308 

Fin3  (P) 

0  0001 

-0.0009 

-0.0054 

Fin4  (L) 

IMi 

-0  0093 

mmm 

Total  fin 

mmm 

0.1203 

Normal  Force 

Finl  (S) 

.  0.3017 

Fin2  (W) 

0.0011 

■I 

Fin3  (P) 

-0.0308 

0.0341 

suited 

0.4132 

Fin4  (L) 

-0.0001 

0.0002 

0.0019 

Total  fin 

0.0000 

0.1550 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward 


x 
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Table  D-3.  Aerodynamic  coefficients  on  planar  fins,  8  =  0°,  Mach  3.0. 


a  0 

2 

4  6 

8 

10 

Axial  Force 

Finl  (S) 

13HM 

■SEEM 

■EMM 

Fin2  (W) 

■tEEM 

MMVE 

bohh 

0.0181 

0.0181 

■EifflESE 

\i\ 1 1  "HM 

0.0158 

0.0150 

0.0152 

EZK9H 

0.0697 

0.0695 

0.0699 

0.0706 

0.0706 

0.0706 

Side  Force 

Finl  (S) 

0.0000 

0.0001 

0.0006 

0.0011 

Eimiw 

Fin2  (W) 

0.0000 

0.0000 

0.0000 

0.0000 

MtHHHIlM 

Fin3  (P) 

0.0000 

-0.0001 

-0.0006 

-0  0011 

-0.0016 

-0.0020 

Fin4  (L) 

0.0000 

0.0000 

-  o.oooi 

Total  fin 

■SEEM! 

0.0000 

0.0000 

■SES53I 

0.0000 

Normal  Force 

Finl  (S) 

0.0000 

0  0382 

0.0776 

0.1179 

0.1551 

0.1881 

Fin2  (W) 

0.0000 

0  0003 

0.0007 

0  0010 

0.0014 

0.0018 

mtmrn 

0.0000 

0.0382 

0.0776 

0.1179 

0.1551 

0.1881 

mssmm 

0.0000 

0.0002 

0  0002 

0  0000 

-0.0004 

-0.0009 

nai 

0.0000 

0.0769 

0.1560 

0.2368 

0.3111 

0.3771 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward 

Table  D-4.  Aerodynamic  coefficients  on  planar  fins,  8  =  10°,  Mach  3.0. 

a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Finl  (S) 

0.0173 

BH 

Fin2  (W) 

0.0180 

Fin3  (P) 

liliiM 

EtlillfiM 

Fin4  (L) 

■nraraai 

0.0148 

0  0150 

Total  fin 

0.0687 

0.0696 

0.0703 

0.0704 

Side  Force 

Finl  (S) 

-0.0001 

0.0006 

0.0011 

mam 

■He! 

Fin2  (W) 

0.0086 

0.0043 

0.0035 

EH 

mmmm 

0.0001 

■TOM 

-0.0015 

bsehI 

-0.0086 

UBiIsFIcEE 

EEH 

Total  fin 

0.0000 

-0.0080 

-0.0082 

0.0038 

Normal  Force 

Finl  (S) 

0.0086 

0.0443 

0.0838 

0.1230 

0.1556 

0.1828 

Fin2  (W) 

0  0001 

0  0004 

0.0010 

0.0014 

Fin3  (P) 

jjfUlfp 

KBi 

Fin4  (L) 

11 

Total  fin 

0.0000 

0.0742 

0.2384 

03119 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 
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Appendix  E.  Force  Coefficients  on  Grid  Fins 


Figure  E-l.  Grid  fin  normal  force,  5  =  0°,  Mach  1.5. 


Figure  E-2.  Grid  fin  normal  force,  5  =  0°,  Mach  3.0. 
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Figure  E-10.  Grid  fin  axial  force,  8  =  0°,  Mach  3.0. 
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Table  E-l.  Aerodynamic  coefficients  on  grid  fins,  8  =  0°,  Mach  1.5. 


a 

0 

2 

4  6 

8 

10 

Axial  Force 

Finl  (S) 

0.0804 

0.0807 

0.0839 

0.0863 

0.0870 

0.0879 

Fin2  (W) 

0.0804 

Mi SMI 

0.0804 

0.0839 

MX'\m 

ESSIII1FSMM 

0.0844 

BEE-iSM 

Total  fin 

■Hi 

0.3498 

0.3524 

Side  Force 

MiggffM 

Fin2  (W) 

-0.0011 

Fin3  (P) 

Bsunval 

-0.0423 

-0.0004 

0.0001 

0.0002 

Total  fin 

-0.0001 

0.0005 

0.0012 

0.0001 

-0.0003 

0.0004 

Normal  Force 

Finl  (S) 

0.0003 

0.0403 

0.0794 

Fin2  (W) 

Fin3  (P) 

Fin4  (L) 

0.0194 

0.0319 

0  0177 

Total  fin 

-0.0001 

0.1097 

0.2069 

0.3061 

0.4174 

0.5218 

S  -  starboard;  W  =  windward;  P  =  port;  L  =  leeward 


Table  E-2.  Aerodynamic  coefficients  on  grid  fins,  8  =  10°,  Mach  1.5. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Finl  (S) 

0  0853 

0.0875 

Fin2  (W) 

0.0821 

BEESIl 

Fin3(P) 

0.0821 

0.0796 

■Hu-231 

Fin4  (L) 

0.0821 

MiHrM 

Total  fin 

0.3284 

■SES 

MEZE3I 

Side  Force 

Finl  (S) 

fcihtira^ 

■HI 

Fin2(W) 

Mil ira 

Fin3  (P) 

-0.0169 

-0.0324 

EjEES 

-0.0321 

-0.0338 

Fin4  (L) 

-0.0117 

-0.0311 

Total  fin 

0.0001 

-0.0267 

EiTiW 

Normal  Force 

Finl  (S) 

MUSI 

WETEgEM 

ItlbltU 

0.2364 

Fin2  (W) 

KftMiciJi 

MMdrM 

0.0714 

Fin3  (P) 

-0.0116 

0.0171 

0.0606 

0.1021 

0.1529 

Fin4  (L) 

0.0170 

0.0413 

Total  fin 

grfnru 

02171 

■55331 

S  =  starboard;  W  -  windward;  P  =  port;  L  =  leeward. 
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Table  E-3.  Aerodynamic  coefficients  on  grid  fins,  8  =  0°,  Mach  3.0. 


a  0 

2 

4 

6  8 

10 

Axial  Force 

OHM 

BVftkWil 

0.0841 

0.0840 

iiUM 

Fin2  (W) 

0.0916 

00939 

— M 

■tliEESl 

0.0841 

W&mEM 

ESS339HH 

fcit/iiiM 

■H 

0.0649 

Total  fin 

0.3217 

0.3228 

0.3263 

0.3269 

0.3247 

0.3183 

Side  Force 

Finl  (S) 

0.0082 

0.0118 

0.0160 

0.0199 

0.0257 

0.0343 

-0.0006 

0 

BMEHM 

-0.0163 

PiTiKM 

Fin4  (L) 

0.0005 

0.0005 

0.0002 

0.0003 

0.0003 

Total  fin 

0.0000 

-0.0002 

-0.0005 

-0.0006 

-0.0007 

-0.0006 

Normal  Force 

Finl  (S) 

0.0385 

0.0795 

0.1206 

0.1598 

01965 

EEE3  I 

0.0247 

0.0783 

wmm 

Fin4  (L) 

0.0082 

0.0203 

0.0214 

0  0104 

-0.0078 

-0.0371 

Total  fin 

0.0000 

0.1049 

0.2060 

0.2962 

0.3755 

0.4353 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 


Table  E-4.  Aerodynamic  coefficients  on  grid  fins,  8  =  10°,  Mach  3.0. 


a 

0 

2 

4  6 

8 

10 

Axial  Force 

Finl  (S) 

0.0838 

0.0798 

0  0834 

0.0840 

0.0821 

WfiEEfli 

Fin2  (W) 

0  0838 

0  0900 

0.0924 

0.0942 

0.0960 

Fin3  (P) 

0.0838 

0.0755 

0.0839 

0.0845 

0  0813 

Fin4  (L) 

0.0838 

0.0713 

0.0652 

0.0636 

0.0630 

0  0622 

Total  fin 

0.3351 

0.3166 

0.3249 

0.3262 

0.3224 

0.3164 

Side  Force 

Finl  (S) 

— ifiTKW 

Fin2  (W) 

-0.0001 

Fin3  (P) 

-0.0052 

-0.0180 

-0.0184 

-0.0205 

-0.0233 

-0.0277 

Fin4  (L) 

-0.0079 

-0.0136 

BHIIITJCM 

Biiinr-KM 

Total  fin 

0.0000 

-0.0169 

BM9 

■iiuaHh 

Normal  Force 

Finl  (S) 

0.0437 

wmm 

0.1920 

Fin2  (W) 

0.0084 

wmm 

0.0781 

Total  fin 

0.0001 

0.2099 

BE51 

wrt^Wii 

0.4452 

S  =  starboard;  W  =  windward;  P  =  port;  L  =  leeward. 


Appendix  F.  Components  of  Aerodynamic  Coefficients 


Figure  F-2.  Components  of  normal  force  for  planar  fin  case,  8  =  10°,  Mach  3.0. 
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Figure  F-3.  Components  of  normal  force  for  grid  fin  case,  8  =  10°,  Mach  1.5. 


Figure  F-4.  Components  of  normal  force  for  grid  fin  case,  8  =  10°,  Mach  3.0. 
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Figure  F-5.  Components  of  side  force  for  planar  fin  case,  8  =  10°,  Mach  1.5. 


Figure  F-6.  Components  of  side  force  for  planar  fin  case,  8  =  10°,  Mach  3.0. 
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Figure  F-7.  Components  of  side  force  for  grid  fin  case,  S  =  10°,  Mach  1.5 


Figure  F-8.  Components  of  side  force  for  grid  fin  case,  8  =  10°,  Mach  3.0. 
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Figure  F-9.  Components  of  axial  force  for  planar  fin  case,  8  =  10°,  Mach  1.5. 


Figure  F-10.  Components  of  axial  force  for  planar  fin  case,  8  =  10°,  Mach  3.0. 
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Table  F-l.  Components  of  aerodynamic  coefficients,  planar  fin  case,  8  =  0°,  Mach  1.5. 


a 

0 

2  4 

6  8 

10 

Axial  Force 

Body 

0.2332 

0.2356 

0.2401 

0.2446 

0.2491 

0.2524 

Total  canard 

0.0245 

0.0246 

0.0247 

0.0246 

0.0244 

0.0241 

Total  fin 

0.0826 

0.0828 

0.0844 

0.0864 

0.0878 

0.0880 

Total 

0.3403 

0.3429 

0.3492 

0.3556 

0.3613 

0.3644 

Side  Force 

Body 

0.0000 

-0.0001 

-0.0002 

-0.0003 

-0  0005 

-0.0006 

Total  canard 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Total  fin 

0.0000 

0.0000 

0.0000 

-0.0001 

-0.0001 

-0.0002 

Total 

0.0000 

-0.0001 

-0.0002 

-0.0004 

-0.0005 

-0.0008 

Normal  Force 

Body 

0  0000 

0.1840 

0.3926 

0  6256 

0.8788 

1.1577 

Total  canard 

0.0000 

0.0350 

0.0689 

0.1013 

0.1316 

0.1598 

Total  fin 

0.0000 

0.1670 

0.3433 

0.5228 

0.6972 

0.8575 

Total 

0.0000 

0.3860 

0.8048 

1.2497 

1.7076 

2.1750 

Table  F-2  Components  of  aerodynamic  coefficients,  planar  fin  case,  5  =  10°,  Mach  1.5 


a  0 

2 

4  6 

8 

10 

Axial  Force 

Body 

0.2452 

Total  canard 

0  0590 

Total  fin 

Total 

w?er~w. 

BMH1M 

Side  Force 

Body 

0.0000 

0.1382 

0  2855 

0.4501 

0.5842 

0.6538 

Total  canard 

0.0000 

-0.0019 

-0.0051 

13MESB 

Total  fin 

0.0000 

-0.0338 

BiTtcgf 

■IWWIKW 

Total 

0.0000 

0.1025 

0.2259 

0.4099 

0.5971 

0.7670 

Normal  Force 

Body 

0.4349 

Total  canard 

0.0599 

— 

Total  fin 

0.0000 

0.3369 

Total 

0.0000 

msnm 
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Table  F-3.  Components  of  aerodynamic  coefficients,  planar  fin  case,  8  =  0°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0.2287 

0.2316 

0.2388 

0.2459 

0.2530 

0.2605 

Total  canard 

K2523I 

■211111 

Total  fin 

0.0699 

0.0706 

Total 

0.3215 

0.3243 

0.3318 

0.3396 

0.3469 

0.3545 

Side  Force 

Body 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0001 

Total  canard 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Total  fin 

0.0000 

0.0000 

0.0000 

0.0000 

0  0000 

0.0000 

Total 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0001 

Normal  Force 

Body 

mm 

0.4979 

Total  canard 

EH 

0.0455 

Total  fin 

i—l 

0  3771 

Total 

1.6401 

Table  F-4.  Components  of  aerodynamic  coefficients,  planar  fin  case,  8  =  10°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0  2309 

0.2337 

0.2404 

0.2473 

0.2538 

0.2601 

Total  canard 

■wig vSm 

llaigB 

Total  fin 

Total 

0.3499 

0.3495 

0.3573 

■lEUjgl 

Side  Force 

Body 

0  0000 

0.0766 

0.1377 

0.1664 

0  1322 

0.0088 

Total  canard 

0.0000 

-0.0096 

Total  fin 

0  0000 

-0.0082 

mm 

■Oukty 

Total 

0.0000 

0.0636 

0.1200 

0.1563 

0.1383 

0.0215 

Normal  Force 

Body 

MB 

0.5048 

0  7798 

1.1997 

0.0320 

0.0481 

0  0644 

0.0807 

Total  fin 

Sf?F7Tai 

0  1561 

0.2384 

0.3119 

0  3770 

Total 

■iXUlIlUM 

0.4936 

0.7912 

1.1561 

1.6574 
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Table  F-5.  Components  of  aerodynamic  coefficients,  grid  fin  case,  8  =  0°,  Mach  1  5 


a  0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0.2372 

0.2393 

0.2440 

0.2485 

0.2525 

0.2554 

Total  canard 

0.0245 

0.0246 

0.0247 

0.0246 

0.0244 

0.0241 

Total  fin 

0.3217 

0.3237 

0  3347 

0.3430 

0  3498 

0  3524 

Total 

0.5834 

0.5875 

0.6034 

0.6161 

0.6267 

0.6319 

Side  Force 

Body 

0.0000 

-0.0012 

-0.0018 

-0.0009 

-0.0012 

-0.0017 

Total  canard 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Total  fin 

-0.0001 

0.0005 

0.0012 

0.0001 

-0.0003 

0.0004 

Total 

-0.0001 

-0.0007 

-0.0006 

-0.0008 

-0.0015 

-0.0013 

Normal  Force 

Body 

0.0000 

0.1347 

0  2896 

0.4695 

0  6773 

0.9239 

Total  canard 

0.0000 

0.0350 

0.0689 

0.1013 

0.1316 

0.1598 

Total  fin 

-0.0001 

0.1097 

0.2069 

0  3061 

0  4174 

0  5218 

Total 

-0.0001 

0.2794 

0.5654 

0.8769 

1.2263 

1.6055 

Table  F-6.  Components  of  aerodynamic  coefficients,  grid  fin  case,  5  =  10°,  Mach  L5. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0.2424 

0  2422 

0.2431 

0.2458 

0.2489 

0.2513 

Total  canard 

0.0637 

0.0634 

0.0623 

0.0608 

0  0590 

0.0574 

Total  fin 

0.3284 

0.3253 

0.3221 

0  3322 

0.3425 

0.3455 

Total 

0.6345 

0.6309 

0.6275 

0.6388 

0.6504 

0.6542 

Side  Force 

Body 

0  0000 

0  1497 

0.3150 

0.4728 

0.5820 

0.5937 

Total  canard 

0.0000 

-0.0019 

-0.0036 

-0  0051 

-0  0062 

-0.0071 

Total  fin 

0.0001 

-0.0267 

-0.0422 

-0.0288 

0.0018 

0  0547 

Total 

0.0001 

0.1211 

0.2692 

0.4388 

0.5775 

0.6413 

Normal  Force 

Body 

0.0000 

0.1528 

0.3406 

0.5281 

0.7390 

0.9879 

Total  canard 

0.0000 

0.0298 

0.0599 

0.0904 

0.1215 

0.1519 

Total  fin 

0.0002 

0.1015 

0.2171 

0.3357 

0.4550 

0.5625 

Total 

0.0002 

0.2841 

0.6176 

0.9542 

1.3156 

1.7023 
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Table  F-7  Components  of  aerodynamic  coefficients,  grid  fin  case,  8  =  0°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0.2341 

0.2371 

0.2441 

0.2514 

0.2581 

0.2647 

Total  canard 

0.0231 

mmm 

— 

Total  fin 

0.3263 

mmm 

■ mm 

Total 

0.5789 

0.5831 

0.5936 

0.6014 

0.6060 

0.6064 

Side  Force 

Body 

0.0000 

0.0004 

0.0005 

0.0004 

0.0000 

Total  canard 

0.0000 

0.0000 

0.0000 

0.0000 

Total  fin 

0.0000 

-0.0006 

-0.0007 

-0.0006 

Total 

0.0000 

0.0002 

-0.0001 

-0.0001 

-0.0003 

-0.0006 

Normal  Force 

nssflBMHI 

0.0000 

0.1356 

0.2890 

0.4846 

Total  canard 

0.0000 

0.0151 

0.0303 

0.0455 

Total  fin 

0.0000 

Total 

0.0000 

fiTEEESM 

Table  F-8.  Components  of  aerodynamic  coefficients,  grid  fin  case,  8  =  10°,  Mach  3.0. 


a 

0 

2 

4 

6 

8 

10 

Axial  Force 

Body 

0.2367 

0.2387 

0.2458 

0.2526 

0.2586 

0.2644 

Total  canard 

0.0470 

0.0471 

0.0473 

wMzSM 

0  0483 

0.0489 

eebhhi 

\mM*m 

ii^a 

03164 

g^ohhi 

0.6297 

Side  Force 

Body 

0  0000 

0.1348 

0.0000 

-0.0095 

-0.0139 

0  0000 

-0.0165 

-0.0031 

0.0256 

0  0511 

Total 

0.0000 

0.0559 

0.1088 

0.1460 

0.1366 

0.0339 

Normal  Force 

Body 

0.0000 

0.7618 

0.0000 

■mb 

0.0644 

Total  fin 

0.0001 

WHi 

UhlfiKM 

| 

Total 

0.0001 

0.8414 

bezmi 
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List  of  Abbreviations  and  Symbols 


A  cell  face  area,  m2 

cal.  caliber  (1  caliber  =  D) 

Q  rolling  moment  coefficient 

Cm  pitching  moment  coefficient 

Cn  yawing  moment  coefficient 

Cp  pressure  coefficient 

Cx  axial  force  coefficient 

Cy  side  force  coefficient 

Cz  normal  force  coefficient 

D  missile  base  diameter,  m 

E  total  energy,  J 

F  inviscid  flux  vector 

G  viscous  flux  vector 

H  vector  of  source  terms 

i,  j,  k  Cartesian  unit  vectors 

M  Mach  number 

MRP  moment  reference  point 

p  pressure,  N/m2 

q  heat  flux  vector 

u,  v,  w  velocity  components  in  x,y,z  directions,  m/s 
V  cell  volume,  m3 

v  velocity  vector  (=  ui  +  vj  +  wk) 

Xq,  location  of  center  of  pressure 

W  vector  of  conservative  variables 

x,  y,  z  axial,  horizontal,  and  vertical  body  axes 
a  angle  of  attack,  degree 

5  canard  deflection  angle,  degree 

v  kinematic  viscosity,  m  /s 

p  density,  kg/m3 

x  viscous  stress  tensor 
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